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Abstract 

We calculate the quantum noise limited displacement sensitivity of a Michelson-Fabry-Perot 
(MFP) with detuned cavities, followed by phase-sensitive homodyne detection. We show that the 
standard quantum limit can be surpassed even with resonant cavities and without any signal- 
recycling mirror nor additional cavities. Indeed, thanks to the homodyne detection, the output 
field quadrature can be chosen in such a way to cancel the effect of input amplitude fluctuations, 
i.e., eliminating the force noise. With detuned cavities, the modified opto-mechanical susceptivity 
allows to reach unlimited sensitivity for large enough (yet finite) optical power. Our expressions 
include mirror losses and cavity delay effect, for a realistic comparison with experiments. Our study 
is particularly devoted to gravitational wave detectors and we consider both an interferometer with 
free-falling mirrors, and a MFP as readout for a massive detector. In the latter case, the sensitivity 
curve of the recently conceived 'DUAL' detector, based on two acoustic modes, is obtained. 
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I. INTRODUCTION 



The sensitivity of interferometers used for the measurement of strain or displacement 
is commonly referred to a so-called standard quantum limit (SQL), calculated considering 
independent fluctuations of the radiation pressure acting on the sensing mirrors and of 
the detected light. Both noise terms are derived from the quantum fluctuations of the 
electromagnetic field. Several studies, starting at least from the beginning of the 80ies, 
have shown that an apparatus can beat the SQL, and accurate quantum calculations are in 

n 

general necessary to find the actual sensitivity. The seminal work by Unruh shows that the 
SQL can be surpassed if quantum correlation characterizes the measuring electromagnetic 
field, and Jaekel and Reynaud ^ show that in this case an ultimate limit is imposed by 
the dissipative part of the mechanical susceptivity. Several schemes have been proposed to 
reach the goal of a sensitivity beyond the SQL, often using additional optical cavities, such 
as the quantum locking 3|, the local meter , 0] , the detuned signal recycling (s 
firstly in Ref. [7| and later analyzed with a deeper attention to quantum noise in Refs. 
It should be noticed that a simple detuning from resonance of an optical cavity allows to 
rotate the field quadratures [lO|, lUl and create a correlation between amplitude and phase 



fluctuations, that are commonly related respectively to radiation pressure and detected field 
fluctuations [12]. This effect is exploited in several proposals of schemes for the generation 



of ponderomotive squeezing 



13 



14 



15| and for quantum non-demolition measurements 



16|. 



A recent work by Arcizet et al. 17|] clearly explains that a detuned Fabry-Perot cavity can 
indeed provide a sensitivity well beyond the SQL, with a frequency behavior very similar to 
that foreseen for interferometers with signal- recycling mirror j^, ^. 

Most of the mentioned studies have been stimulated by the development of large interfer- 
ometric detectors of gravitational waves (gw). Recently, a new class of gw detectors has been 
conceived, based on huge masses kept at cryogenic temperature and called DUAL detectors 

3, 



IS 



20l |. Differently from previous massive cryogenic antennas (such as Weber bars). 



the DUAL system do not exploit particularly a mechanical resonance of a solid body, but 
it takes advantage of elastic forces to achieve a useful sensitivity in a wide frequency range. 
At this purpose, it has to give up to the usual resonant mechanical amplifier, and it needs 
a very sensitive readout. One possibility is using a Michelson interferometer with suitable 
Fabry-Perot cavities in the two arms (Michelson- Fabry- Perot, MFP) 2l|. The readout would 
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be in principle similar to the large gw interferometers like VIRGO 22| and LIGO 23| , but a 
more complicated mechanical susceptivity and response function to gw must be taken into 
account. The definition of SQL is less obvious than in usual interferometric detectors, as 
well as the possibility to surpass it, and a study of such a system fully including quantum 
noise is still lacking. 

In this article we calculate the sensitivity of a MFP interferometer with detuned cavities. 



The calculation is very similar to the one described in Ref. [17| for a simple cavity, and we 
find indeed spectral curves well beyond the SQL, with shapes comparable to those typical of 
signal-recycled interferometers. In addition, we include cavity losses for a better comparison 
with realistic experimental schemes, we introduce as additional degree of freedom the choice 
of the final detected field quadrature, and we apply the results to both a standard free-falling 
masses interferometer and to a DUAL detector. 



II. THEORETICAL MODEL 



A. Simple cavity 



Before describing our complete model, we analyze the paradigmatic case of a Fabry- 
Perot cavity with a movable mirror, neglecting mirror losses and cavity field delay (short 
cavity regime). Such a calculation is reported in details in Ref. l?!, and we only add the 
choice of the detected field quadrature that can be performed by using a local oscillator with 
tunable phase. Such simplified scheme is useful to understand the physical meaning of the 
phenomena that will be observed in the complete system. 

We use the semi-classical formalism described in Ref. [l^, valid in the limit of strong 
fields, where quantum field fluctuations are treated as classic stochastic variables. 

In the limit of high finesse and nearly resonant conditions, the equation for the cavity 
field a reads 

(-7 + #)a + ^am = (1) 

where if) = 2kL Mod[27r] is the phase detuning from the closest resonance {k is the laser field 
wavenumber and L is the cavity length), 27 is the input mirror intensity transmission, 
is the input field. The electric fields are normalized such that |ap is a flux of photons. 
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The input /output coupler boundary conditions are 

Oiout = -o-in + a/27 a 



(2) 



In a linearized analysis, the general electric field can be considered as a sum of a steady 
state a (which is null for vacuum fields) and the fiuctuations 5a{t) around it. The steady 
state of the intracavity and refiected fields are respectively 



and 



- _ 

a = rj otin (3j 

— lip 



(Xout = 7j CXin ■ (4j 

— lip 



The steady state of the cavity detuning ip is 

-i/i = 7/;o + 4M\|ap. (5) 

Here ipQ is the cold-cavity detuning (for vanishing laser field) and the last term in Eq. ([5]) is 
the radiation pressure effect, where x is the movable mirror susceptivity. 

The linearized fiuctuations 5il){t) of ip around its steady state and the field fiuctuations 
can be written in the Fourier space defining 5a{t) = 6a{Q)e~^^^ , 6a* (t) = (5a*(fi)e~*™ and 
Sip{t) = 6ip{Q)e'^^^. The equations for such fiuctuations read 

d-pin) = 5^0 + 4 hk\{a*6a + a6a*) (6) 

(7 — iilj)Sa — iaS'ip = \/2rfSain (7) 

where 6ipo is the signal to be detected, and the equation for 6a* is the conjugate of Eq. ([7j). 
It is useful to use the quadratures of the field fiuctuations, defined as 

6p = 6a + 6a* ; 6q = i{6a* — 6a) . (8) 

6p and 6q correspond respectively to the amplitude and phase fiuctuations, referred to the 
input mean field that is taken as real (i.e., q;*„ = ain)- 
Using such quadratures, the equation (JTj) becomes 



'J 6p + ip 6q — i{a — a*)6ilj = 6pin (9) 
'y 6q - ■p 6p - {a + a*)6'ijj = ^/2rf 6qin ■ (10) 



To write clearer expressions, we define \l/ = 7/^/7 (detuning normalized to the half cavity 
linewidth) and we use a normalized input laser power p (with the dimensions of a force 
divided by a length, i.e., the inverse of a susceptivity) defined by 

p=16 -^ = 16^ 11 

7^ 7^c 

where Pm = hkcal^ is the real input power (c is the speed of light). With these definitions, 
the expressions ([3]) and (jll) for the fields steady states become 

tt=\/ 7=== "in (12) 

aout = e^'"^ oiin (13) 

where 77 = arctan\I'. 

The quadrature fluctuations, according to Eqs. ([9]) and ffTOj) . are rotated by the same 
angle as the average field. Indeed, we have inside the cavity 



Sp\ [21 ( cos 1] —smr]\ I 6pin \ / 2 \ ^ I —sm2r] 




1 sin T] cos T] I \ hvn J V7/ 1 + COS 277 



5i) (14) 



where we have used Eq. (IT^ for replacing a. 

For the reflected fields, we find relations similar to Eq. ( fT4l) . where the rotation angle for 
the field quadratures is again the same as for the steady state field: 

cos2r7 -sin2r/\ / 5pin \ ^ 4 a^^. / -sin2r7 ^^^^ 
sin 2?7 cos 2r] I \ 6qin I 7 1 + ^1/^ \ cos 27] 

The equation for Sip, in the simplified notation, reads 

Sijj = d-i/jo + 4 /i/e^x-^/" y==|[cosr7 5p + sinr] 5g] (16) 





H 







and, replacing Eq. (fT4|) for (5p and 5g, we obtain 



2 ^ _ ~ 

^Pi^ 1 I ,T,2 "^'^ 

7 1 + 



(17) 



-71 + ^2 

The detected quadrature can be chosen at will by tuning the homodyne angle w, and the 
corresponding fluctuations are 

6Eout = Spout cos w + 6qout sin w . (18) 
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We remark that, with respect to the input field, an amphtude detection is obtained for 
w = while a phase detection corresponds to w = 7r/2. On the other hand, with respect 
to the output field (see Eq. (IT^ ). w = 2rj corresponds to a pure amplitude detection and 
w = 21] + 7i/2 to a pure phase detection. With w = 2ri we have SEout = ^Pm while for 
w = 21] + tt/2 we obtain 



4 a 
7 1 + ^' 

In general, for a given homodyne phase w, the detected field fiuctuations are 



4 a ■ 

6Eout = Spin cos{w -2r])+ [ 6qin + - ^ 6ip ] sm{w - 2t]) . (20) 



Replacing the expression ( |T71) for 6ip in Eq. ( l20l) . we get 

SEout = Spin cos(w - 2r]) + sm{w - 2i]) 
with 



4 6ipo 2 A 



(21) 



-4= n , .,„,, PX • (22) 



The last term in Eq. fl^Tl) is the contribution of the radiation pressure noise, proportional 
to the input field amplitude fiuctuations 6pin. For w = 21] + tt/2, the amplitude fiuctuations 
only enter into the detection noise through such radiation pressure term. Changing the 
detection angle from {w = 21] + tt/2), the contribution of the phase fiuctuations 6qin to 
the signal-to-noise ratio (SNR) remains the same (both 6qin and 6ipo multiply the angle- 
dependent factor sin(w — 2i]) ). On the other hand, a further contribution of the amplitude 
fiuctuations comes into play, that can compensate the radiation pressure fiuctuations. The 
optimal SNR is obtained when the contribution of 6pin is completely cancelled, a situation 
occurring for an angle Wopt such that 

, . 1 + A . . 

tan{Wopt - 2i]) = — — . (23) 

In this case, the SNR is only limited by the phase fiuctuations, and it increases with the 
coefficient of tpQ within square brackets in Eq. fl2T|) . In particular, the sensitivity increases 
for negative detuning (i.e., for A < 0) and one can get in principle unlimited SNR if the 
denominator of this coefficient vanishes, i.e., for A —1. This can happen for large enough 
power and/or susceptivity, more precisely if px > 16/3^/3 . Such increase in sensitivity at 
negative detuning is usually interpreted as due to a modified effective susceptivity ('optical 
spring') originated by the position-dependent radiation pressure force. 




FIG. 1: Left: the optical configuration considered is a Michelson- Fabry- Perot with the additional 
free choice of the detected field quadrature. Right: mechanical scheme, with self- and cross- 
susceptivities. 

B. Complete system 

The optical scheme we consider, reported in Fig. ([1]), is a MFP with the addition of an 
homodyne balanced detection at the output. In this way, one can choose the quadrature 
of the output field to be detected as in the paradigmatic case just analyzed. Vacuum 
fluctuations are introduced through the mirrors losses in the cavities and through the output 
port of the Michelson beam splitter, while laser field fluctuations enter through the beam 
splitter input port. 

Concerning the mechanics of the system, we consider a susceptivity matrix to include 
the possibility of changing the length of one arm by acting with a force on the other arm. 
Such a possibility is important in the case of an interferometer mounted on a solid body. A 
sketch of the mechanical scheme is included in Fig. ([T]) (right). 

The equations for the electric fields in the two cavities (labeled by i,j = (1, 2)) read 

T-i^ = -(% - + V'2^af + ^/2^i6aH (24) 

where = 2Lj/c is the cavity roundtrip time, c is the speed of light, 27^1 is the input 
mirror intensity transmission, 27;, are the roundtrip intensity losses (including transmission 
from the back mirror, absorption and scattering in both mirrors, diffraction losses, etc.), 
li = Imi + In, Sc(ii are the vacuum fluctuations entering through cavity losses which are 

7 



mimic by a partially transmitting output mirror. 

Assuming an ideal 50% beam splitter, the input fields a*" of the two cavities are 

«i = (25) 

«r = ^^^^^ (26) 

where am is the laser input field and 5av are the vacuum fluctuations entering through the 
beam splitter output port. In Eqs. (12511261) . to simplify the notation, we have neglected the 
phase difference between a™ and introduced by the length difference in the paths from 
the beam splitter to the two input cavity mirrors. This phase difference will be re-considered 
at the output port of the beam splitter (Eq. (l28l) ). 

The cavity input /output coupler boundary conditions are 



-af + V^ai. (27) 



At the output port of the beam splitter, fields are recombined giving: 

aB5 = ^(-«r* + «rv^) (28) 

where the phase 6 accounts for the double path difference between the beam splitter and 
the two input mirrors. Finally, the observed quadrature of the output field can be chosen 
by changing the detection phase w. 

The steady state of the intracavity field is obtained by zeroing the time derivative in 
Eq. ( !24l) . using Eqs. ( !25ll26l) and neglecting the field fluctuations: 

oti = — — am • (29) 

7i - m 

Using Eq. (1271) . we find for the steady state of the reflected field: 

Imi + 111 - ii'i V2 

The cavity length is sensitive to several kinds of forces acting on the system, including 
classic deterministic (e.g., the gravitational wave effect), stochastic (thermal noise) and 
quantum forces (the radiation pressure acting on the mirrors). The steady state cavity 
detuning ipi can be written as 

iJi = ^jjoi + Ah^xllcxil^ + Ahk^llajl^ (31) 



where x^j is the stationary (zero-frequency) susceptivity matrix. 

The equation for the fluctuations of the cavity phase detuning 6ipi{fl) (expressed in the 
Fourier space) is 

6ipi{n) = dijQi + Ahk'^[xii{9){a*6ai + aida*) + Xij{^){a*6aj + ajda])] (32) 

where dipQ contains the effects of thermal and external noise, and of the gw signal and Xij{^) 
is the susceptivity matrix. 

The complete set of equations for the fields and cavities displacements are reported in 
Appendix A. The complete expressions with different parameters for the two cavities are 
useful for a future numerical analysis of the effects of the asymmetries and of the allowed 
tolerances in the parameters. However in this work, for a simpler understanding of the 
physical phenomena, we will take identical cavities with 7^1 = 7m2 = 7m) In = I12 = li, 

ri = T2 = T. 

As a further restriction to our analysis, we will consider that: 

a) the Michelson interferometer working point gives a dark fringe at the beam splitter 
output port. This corresponds to setting equal distances between the beam splitter and the 
two cavities, i.e., e*^ = 1; 

b) the two cavities have the same detuning from the laser frequency: ipi = ip2 = i^] 

c) the mechanical system is symmetric: Xii = X22 = Xs and X12 = X21 = Xc- 

These three conditions determine the cancellation of the effect of the input laser field 
fluctuations in the output. This requirement is important for a system working in the 
acoustic frequency range. Indeed, while it is very difficult to reduce the laser field amplitude 
fluctuations at the shot noise 

limit is even more difficult and far from being demonstrated in strong power laser fields. 

The field fluctuations SEout (seen by the homodyne detection) are described by a vector 
of coefficients Vout, multiplying the input fluctuations Xj„: 

SEoum = Vout ■ X,„(fi) (33) 

with Xin and Vout given respectively in equations and (jHS]) of the Appendix, where we 
have used for convenience the quadratures of the field fluctuations. 

As already remarked, the coefficients V^tit[(5pj„] and that multiply the input 

laser field fluctuations 6pin and 6qin are null in the completely symmetric case that we are 
considering. 
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Besides the previous defined normalized detuning \I/ = ?/'/7, we use for a more compact 
notation Tm = 1m/ 1 (in tlie case of loss- less cavities, Tm = 1); ^cav = i/t' (cutoff angular 
frequency of the cavity) and the normalized input laser power p is now 

P = 16 5 . (34) 

The expression of the sensitivity Sl{Q) in the detection of S{Li — L2) (defined as the signal 
spectral power with unitary signal-to- noise spectral density) is given in Appendix A, Eq. fl^ . 
We see that in SLi^t) the susceptivities only appear as a difference Xs — Xc- We define in 
the following Xs-Xc = X- 



III. GENERAL DISCUSSION 



The first situation that we consider is with the cavities at resonance {ip = 0) and pure 
phase quadrature detection {w = 7r/2). Such configuration corresponds to the present 
operative gw interferometers VIRGO and LIGO, and we will define it in the following as 
'normal case'. The Sl{Q) of Eq. fl^ becomes: 

2" 



Sr 



n 1 + 



n \2 



^cav ' 



V 



v\x\ 



1 + 



^ \2 



(35) 



It can be seen that Si^ can be written as a sum of a 'displacement noise' term S^x and a 
term proportional to a 'force noise' Sff'- 



Sl — Sxx + \x\ Sff 



with 



Sxx S 



xx^FF 



h^Tr 



(36) 



(37) 



The origin of the two terms can be found respectively in the intracavity field phase noise 
and amplitude noise. The first term limits the detection sensitivity of a phase signal created 
by the mirrors displacement; the second term is due to the fluctuations of the radiation 
pressure acting on the cavity mirrors. For each detection frequency Q, the optimal Sl{Q) 
is reached for 5*^.^. = h \x\ / y/T^- With this condition, Sl is equal to what we call standard 
quantum limit: 



SQL{n) 



2n\xm 



(38) 
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With long cavities (i.e., when Vt/Vtcav is not neghgible) and/or a frequency- dependent suscep- 
tivity, the SQL defines an envelop of possible spectral density curves, each one determined 
by the choice of the input power. 

Already with resonant cavities, tuning the homodyne phase w allows to change the sen- 
sitivity significantly. still assumes a simple form: 



+ 




(39) 



cosw . 

sm w 1 + ' -^1^ 



that, for loss-less cavities (F^ = 1) and real x, can be further simplified to 

2' 



i + (^)2 

P 



PX 



cosw 



sin w 1 + (o^)^ 



(40) 



The physical interpretation of Eq. (1401) is simple, as already explained for a simple cavity. 
Like in the 'normal case', for resonant cavities there is no quadrature rotation, the vacuum 
fiuctuations 6pv and 6qv are transferred respectively to the amplitude {6pBs) and phase 
i^lBs) fiuctuations of the output field, and the radiation pressure is proportional simply to 
6pv- The cavity length fluctuations, besides the signal, contain a term due to the radiation 
pressure, and a term proportional to the length fluctuations is present in the output fleld 
phase. As a consequence, the amplitude fluctuations 6pv are found both in the output fleld 
amplitude quadrature SpBs (giving the flrst term in round brackets of Eq. PUI) ) and in the 
output phase 6qBs (second term in round brackets of Eq. (HOl) ). A suitable choice of the 
homodyne phase, and therefore of the detected quadrature, brings these terms to cancel 
each other. Also for 7^ 1, choosing 



w = Wq = arctan 



y/TmpRex 

the contribution of amplitude fluctuations is minimized and the Sl becomes 



(41) 



Sl = SI = S^, + Sff[{1 - rm)Rex^ + Irnx^] ■ 



(42) 



For loss- less cavities (F^ = 1) and real the radiation pressure noise can be completely 
cancelled and we have = Sxx- 
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FIG. 2: Sensitivity 5^ as a function of input power p, for ^cav oo, Tm = 1, and constant, 
real x- Solid line: 'normal case' = 0, w = 7r/2); dashed line: general resonant case = 0, 
w = arctan(— 1/px)); dotted lines: ^ = —0.4 and w = tt/2 (a), w = 2.0 (b), w = 2.3 (c). 

For frequency- dependent susceptivity and/or long cavities, defines the locus of the 
minima of spectral curves that can be tuned by changing w. The value of w necessary to 
achieve the best sensitivity at a chosen frequency is given by Eq. (14T|) . We remark that in the 
'normal case' the SQL is an absolute limit, while in the general resonant case (in particular, 
for F^ = 1) the sensitivity limit given by S^x can be decreased at will by increasing the laser 
power. 

For detuned cavities, both the radiation pressure and the output field fluctuations contain 
a mixture of 6pv and 6qv and the force noise is not due any more solely to the pure amplitude 
fluctuations 6pv- A clear discussion of the phenomena involved in this case is given in 
Ref. [171 . 

For a better understanding of the physics, we analyze the case of a constant, real, positive 
X and very short cavities {Qcav ^ ^) with negligible losses (F^ — 1). The 'normal case' 
sensitivity is shown in Fig. ([2]) with a solid line, as a function of the input power. The 
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SQL is reached for p = psql = ^^t p > Psql the sensitivity is worse due to strong 

radiation pressure effect, for p < psql it is deteriorated by phase noise. The sensitivity 

for the general resonant case is shown in Fig. ([2]) with a dashed hne. In this case, 
coincides with 5*^.^.. For low power it approaches the 'normal case' sensitivity (that is here 
dominated by S^x), but for p > psql/'^ it surpasses the SQL. If we allow for different values 
of the detuning \E', Sl can decrease well below the SQL and even below 5'°, as shown in 
Fig. ([2]) with dotted lines, and it is unlimited if p is strong enough. The physics behind this 
effect was previously explained in Section III A[ In short terms, it is the result of: a) the 
cancellation of the input amplitude fluctuations (thanks to a good choice of the homodyne 
detection angle); b) the modified effective susceptivity, that increases the sensitivity to the 
mirror motion. 

Considering now a complex susceptivity Xi taking into account mechanical dissipation, it 
can be shown that the minimal sensitivity is h Imx, as already found in Ref. ^ for a MFP 
with resonant cavities and squeezed input fields, and by Arcizet et al. [17] for a detuned 
Fabry-Perot cavity. This phenomenon can be understood as follows. In the detected field, 
the fluctuations 5pv and 5qv are present for two different reasons: directly in the field 
reflected by the interferometer (purely optical effect), and because of the length fluctuations 
induced by radiation pressure (opto-mechanical effect). As we have already seen, the same 
situation is found in the case of resonant cavities, but only for 5pv- The two effects give 
different linear combinations of 5pv and with real coefficients depending on ^ and 
w (we are considering Vtcav — ^ oo). However, the opto-mechanical effect is mediated by 
the susceptivity x- If X is real, an appropriate choice of ^ , w and p (and therefore of 
the coefficients multiplying Spy and 6qv) can bring to a complete cancellation between 
the purely optical and opto-mechanical effects. As already explained, the same situation, 
occur ing in the case of resonant cavities only for 6pv, defines S^. However, if x has an 
imaginary component, the radiation pressure fluctuations cannot be completely cancelled. 
In other words, the total detected fluctuations cannot be completely deleted because of the 
de-phasing between intracavity intensity changes and cavity length changes introduced by 
the complex x- In spite of this interesting physical result, a region with significant imaginary 
part of the susceptivity is of limited practical interest and in the following we will assume a 
real x- 
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FIG. 3: Sensitivity Sl in the case of free-falling mirrors, as a function of the frequency 0,, for 
^cav OO, Tm = I, p = l- Solid line: 'normal case' {"^ = 0, w = it/2); dashed line: resonant case 
limit 5^; dash-dotted line: SQL; dotted lines: = and w = 1.2 (a), w = 0.8 (b), w = 0.4 (c), 
w = 0.2 (d). 

IV. FREE-FALLING MIRRORS 

We apply now our results to find the sensitivity of an interferometer with free-falling 
mirrors, still with loss-less short cavities. Such a scheme is a good approximation for mirrors 
suspended to a pendulum with low oscillation frequency. The susceptivity can be written as 

X^-.„(|)^ (43) 

where Xo and Qq are constants. We remark that a finite zero-frequency susceptivity (and 
therefore a defined pendulum oscillation frequency) is necessary in the calculation of the 
steady state and stability of the system, but it can be neglected in the evaluation of the 
spectra. 

As well known, in the 'normal case' for each particular choice of the power Pj„ the SQL 
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is reached at a corresponding frequency ^sql given by 

^Iql = PXo^l , (44) 

where Sl is: 5*^(^5'^^) = 2h/ (y/TmP)- At high frequencies, 5*^ tends to the asymptotic 
value Sl{Q oo) = jSl{^sql), hmited by phase noise, while below ^sql the Sl increases 
as due to radiation pressure fluctuations. An example is shown in Fig. [3] (solid line). 

Keeping resonant cavities but changing w, the sensitivity curves can surpass the SQL, 
with minima lying on the horizontal line given by Sxx (Fig- E]). We remark that in our 
scheme (with completely symmetric cavities), the average field at the output port vanishes, 
therefore the local oscillator power can be kept low. At first order, the balanced homodyne 
detection scheme is not sensitive to the local oscillator noise. Therefore, the phase tuning 
element (e.g., a phase mirror or an electro-optic modulator) is not critical for the noise 
budget. As a consequence, this sensitivity tuning technique can be easier to be implemented 
than other schemes (e.g., tuning the MFP cavities or a signal recycling mirror). 

If we now allow for detuned cavities we see that even S^ can be largely surpassed. Some 
examples are shown in Fig. (jlj) for a detection phase w kept at tt/2. We remark that the 
position of the minimum Sl shifts toward low frequencies at increasing detuning. Therefore, 
the input power should be increased to keep the same optimal frequency range. Similar 



results are described in Ref. [17|] (see their Fig. 3)[25|. 

Tuning the phase w allows to obtain families of sensitivity curves, one for each fixed value 
of detuning, whose envelopes are very broad and low. Some examples are shown in Fig. (Q, 
where the gray solid curve represents the absolute sensitivity limit that can be reached 
optimizing, for each frequency, both the phase w and the detuning. In the inset of the figure 
is shown the optimal phase w, as a function of the frequency, that is necessary to obtain 
the envelope curves. For resonant cavities, inserting the susceptivity (H3l) in equation (HI]), 
we find that the optimal phase is w = Wq = arctan(fi/i75Qi)^. The other curves differ 
from Wq for a very fiat additional phase. We notice that, besides choosing a particular 
detection phase to optimize the detection at the preferred frequency, one can also obtain a 
sensitivity curve corresponding to one of the above envelopes by adding a suitable frequency- 



dependent quadrature rotation at the output 
so-called variational-output interferometers ( [26 



3efore the homodyne detection, like in the 

3)- 



Moderate cavity losses do not change qualitatively the above features, and also quanti- 
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FIG. 4: Sensitivity Sl in the case of free-falling mirrors, as a function of the frequency 0,, for 
^cav OO, Tm = I, p = l- Solid line: 'normal case' {"^ = 0, w = it/2); dashed line: resonant case 
limit Sl; dash-dotted line: SQL; dotted lines: w = Tr/2 and ^ = 0.4 (a), ^' = 0.6 (b), ^ = 0.8 (c), 
^ = 0.9 (d), ^' = 1.1 (e), ^ = 1.2 (f), ^ = 1.4 (g). 

tatively the enhanced sensitivity is well preserved. An example is shown in Fig. ([6]). Even 
at Tm = 0.5 (input mirror transmission equal to the other roundtrip losses) the dip in the 
Sl is well pronounced, below the SQL. 

The expression of 5*^ including the cavity response is reported in the Appendix A 
(Eq. (l65|l ). The simplified expression for the 'normal case' is 

Sl{Q) = SLiQsQL) ■ \ 

Two examples are shown in Fig. ([7]), for Vtcav = O.Sfi^QL and VLcav = '^^sql (dark, thin 
lines). In any case, the SQL remains the same (as in Eq. fl55]) ). and at high frequency (above 
VLcav and ^sql) the S^ now increases as VL^. 

For detuned cavities, the most interesting situation is when VL^av is around or above ^sql- 
In this case, for VL < Qsql the Sl is not significantly modified with respect to the previously 
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^cav 



(45) 




FIG. 5: Sensitivity Sl in the case of free-falling mirrors, as a function of the frequency il, for 
^cav CO, Tm = 1, p = 1. Black solid line: 'normal case' = 0, w = tt/2); long-dashed straight 
line: resonant case limit S^; dash-dotted line: SQL; dashed lines: the phase w is optimized for each 
frequency, ^' = 0.06 (a), ^' = 0.1 (b), ^' = 0.2 (c), ^' = 0.3 (d), ^' = 0.4 (e); gray solid line: both 
w and ^ are optimized for each frequency. In the inset: phase w optimized for each frequency, for 
the values of ^ used in the main plot. 

studied situation of r ~ 0. At high frequency (well above ^cav) we can write 

Sl increases as Q'^ like in the 'normal case', and is always above it, with an optimal detec- 
tion phase of w = I + arctan(\[^). What is more peculiar is the behavior at intermediate 
frequencies. Here a further increase in sensitivity is found (see the 'bump' in the light (read 
online) dotted curve of Fig. ([7])). The Sl is here lower than in the case of negligible r, yet 
remaining above the SQL. 
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FIG. 6: Sensitivity Sl in the case of free-falling mirrors, as a function of the frequency 0,, for 
^cav oo, p = 1. Solid lines: 'normal case' {"^ = 0, w = 7r/2, from lower to upper curve 
Tm = 1, 0.8, 0.5); dashed line: resonant case limit for = 0; dash-dotted line: SQL for 
Tm = 0; dotted lines: w = tt/2, ^ = 0.6, and from lower to upper curve F^ = 1, 0.8, 0.5. 

V. DUAL DETECTOR 

A DUAL gw detector exploits two oscillation modes of a mechanical system with a readout 
symmetric with respect to the center of mass. Do to the geometry, the responses to the 
readout force of the two modes must be summed, while the responses to the tidal force of 
the gw are subtracted. In the frequency region between the two resonance frequencies, the 
susceptivities of the two modes are in anti-phase, giving a reduced response to the readout 
force and an enhanced response to the gw. Different configurations have been proposed 



and studied: two nested spheres 181], where the relevant modes are the first quadrupolar 



mode of the inner and outer sphere; two nested cylinders 19|, again acting on the first 
quadrupolar mode of the nested bodies; a single hollow cylinder, exploiting its first and 
second quadrupolar modes; a symmetric set of cylinders, where the first DUAL mode is 
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FIG. 7: Sensitivity Sl in the case of free-falling mirrors, as a function of the frequency for 
Tjn = 1) P = 1- Thin, dark lines (black online): 'normal case' {"^ = 0, w = 7r/2); thick, light lines 
(red online): ^' = 0.6, w = 1.3. Solid lines: flcav oo; dotted lines: i^cav = '^^SQL', dashed lines: 
^cav = O.Si^sQL- Dash-dotted line: SQL. 



2Q|. 



given by the link between them and the second one by their first oscillation mode 

For our analysis, we consider the simple case of very low (vanishing) first resonance 
frequency. The susceptivity (i.e., the difference between Xs and Xc) can be written as 



1 1 

-7^ + 



Xo^^o (47) 



fi2 ' ^{ni-n^)_ 

where is the (second) DUAL resonance frequency, /i is an adimensional modal mass and 
the frequency is normalized to i^o- A typical choice of Qq is Qq = Vs/R, where is a typical 
dimension of the detector and Vg is the sound velocity in the detector material. Xo depends 
on the material properties and geometrical configuration of the detector. 
The response to the gw is 
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FIG. 8: Sensitivity to a gravitational wave Saw of a DUAL detector, as a function of the frequency 
J7, in arbitrary units. Dark, solid line (black online): SQLqw- Light solid line (red online): 'normal 
case' (resonant cavities, w = it/2) with p = 0.4. Dotted lines: 'normal case' with p = 2 (a), p = 0.1 
(b). 

where ^gw is an adimensional gw sensitive modal mass. Both /i and ^cw are normalized to 
the detector physical mass and can be calculated from the detector geometry: yU is a surface 
overlap integral between the radiation pressure of the field interrogating the surface and the 
considered mechanical mode; iiqw is a volume overlap integral between the quadrupolar gw 
force and the mechanical mode. We take for example in this article the 'QUAD' detector 

n 

described in Ref. [20|, with four parallel molybdenum cylinders. In this case, we have 
/i = 0.28, liGW = 1.1, = 1.4. 

Being interested in the sensitivity to a gw signal, we can define it as 
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In the case of the free-falhng masses MFP considered above, we have Hqw = 1 and therefore 
Sow = Sl- The standard quantum hmit is defined similarly to Eq. fl38|) : 



SQLgw = (50) 

V^r\n|-?/GVK(^)P 

and is shown in Fig. ([8]) with a dark, solid line. 

In the 'normal case', varying the input power gives a family of sensitivity curves, limited 
by the SQLgw- A proper choice of the power gives a rather flat curve. For lower power, 
the sensitivity is peaked around Qr where both the gw signal and the radiation pressure are 
amplified. A deep in the radiation pressure effect is found around a particular frequency Qrp 
where the interference between the two detector modes gives xi^rp) = (for the parameters 
here considered it happens at Qrp — 0.655). For high input power, around this frequency 
we have the best sensitivity (Fig. [HI). 

Also for a DUAL detector, keeping the cavities at resonance and varying w one obtains a 
set of curves whose envelop, shown in Fig. ([9]) with a dashed line, is well below the SQLqw 
and allows to widen the detector sensitivity around Vtji. Changing the detuning \E', we find a 
Sgw getting even lower. As shown in Fig. ([9]), depending on the parameters, the reduction 
below SQLgw can occur either in a range between VL^p and VLr (where x > 0), or below VL^p 
and above VLr (where x < 0). It can be shown as a general property that, for fixed detuning 
and phase w, it is not possible to obtain a sensitivity curve falling below the SQL both in 
frequency regions with positive x and in regions with negative x- 



VI. CONCLUSIONS 



We have presented the analysis of a Michelson-Fabry-Perot interferometer with the addi- 
tion of the free choice of the detected field quadrature. Our study fully accounts for quantum 
noise, including radiation pressure effects and possible losses in the cavities, and we consider 
in particular the case of detuned optical cavities. 

The use of a susceptivity matrix allows to extend the applicability of our results from 
the usual, free-falling mirrors to the most general mechanical system, including any kind of 
DUAL detectors with mirrors installed on elastic, mechanically coupled test masses. In view 
of the development of extended numerical evaluations, considering construction tolerances 
and variations of system parameters, we have given complete general expressions of the 
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FIG. 9: Sensitivity to a gravitational wave Sqw of a DUAL detector, as a function of the frequency 
$7, in arbitrary units. Dark, solid line (black online): SQLqw: the other curves are for p = 0.4; 
dashed line: resonant case limit; light solid line (red online): 'normal case'; dotted lines: w = tt/2 
and ^' = -0.5 (a), ^' = 0.5 (b). 

expected output signals. 

A physical analysis of simplified expressions shows that, thanks to the possible choice of 
the detected quadrature, the sensitivity can surpass the SQL even with resonant cavities. 
This is due to the complete cancellation of the radiation pressure fluctuations in the detected 
output field quadrature. With detuned cavities, in the case of real susceptivity and loss-less 
cavities, the peak spectral sensitivity is unlimited if the susceptivity and/or the laser power 
are strong enough: the ultimate sensitivity is only limited by the reactive component of 
the susceptivity. We also show that reasonable optical losses do not critically modify the 
performance. 

The scheme that we consider, with the choice of the detected field quadrature, is easier to 
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be implemented than other proposed and tested configurations (e.g., the resonant sideband 



extraction 



28l | or the signal-recychng cavity 



It must be remarked, however, that here 



the sensitivity enhancement occurs in a region of high susceptivity (corresponding in the 
case of free-faUing mirrors to the low-frequency region). 

While, for a chosen value of the detection phase, the sensitivity is optimized at a particular 
frequency (i.e., for a particular value of the susceptivity), one can implement a frequency- 
dependent rotation of the output field before detection, as proposed and analyzed in several 
works [26|, |27|, in order to enlarge the useful bandwidth. Without discussing and detailing 



the possible experimental schemes suitable at this purpose, we give the obtainable sensitivity 
curves that are indeed broad and deep with respect to the SQL. 

The use of detuned cavities is again somehow simpler than other schemes that are based 
on additional mirrors and cavities or on the use of squeezed light, yet allowing to observe 
interesting noise reduction effects similar to those reported in the literature concerning the 
mentioned configurations. We remark that detuned cavities imply in general critical stability 
problems due, e.g., to radiation pressure and photo-thermal effectsQ, The analysis of 
this topic is beyond the purpose of this article, and it should be developed for each particular 
system including the full frequency response and the possible active stabilization. 



VII. APPENDIX A 



The equations for the field and cavity phase linearized fiuctuations in the Fourier space 
can be written in a compact matrix form. Defining an input fiuctuations vector y4j„(fi) and 
a cavity fiuctuations vector A{Q) as 



A{n) 



^m(^) 



Sam 

equations (12411261) . fl32|) can be written as 



\ 5lp02 ) 



(51) 
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(52) 



where 



M{Q) 



( —itVL — iipi + 7i 





V 4 /c^ X21 ai 



—iai 

-i r i7 + i-i/^i + 7i 2 a* 

Ak'^hxiiOi 1 4k'^hxi2C(2 

— 2 r — iil!2 + 72 



4 7i ;>^2i at Ak'^hx22a*2 






Ak'^hXu "2 




-i r n + i-?/'2 + 72 

4 /C^ X22 0.2 

(53) 

It is useful to write the field fiuctuations in terms of amplitude and phase quadratures, 
defined as 

5p = (5a + bo* ; 6q = i{6a* — 6a) . 



\ 




-ia2 
i 



Eq. (152|) becomes 

N(fi) ■ x{n) = xu^) 

with fiuctuation vectors X{Q) and Xj„(f2) and system matrix N(ri) given by 



(54) 



(55) 



x{n) 



\S^2j 



V^^Pl2 + i^Pin + ^Pv) 

^'4'02 



(56) 



-iTiil + 7i 






-iTiil + 7i 






2ain^Jml4'l 













-4fe^^«i„7i^7mlX2l -4A;^a«i„^7mlX2i'/'i 



7l2+^ 



1 








-4/c^aai„72v'7m2Xl2 -4fc^ V')'m2 Xl2 V'2 



72^+V'2 



-2X2^ + 72 
-^2 



72^+'/'i 

^2 

-2X2 + 72 



-4fc^7),ai„72 V7m2X2 -4fc^S,ai„ ^7m2X2'/'2 



72^+V^ 



72^+^ 



■(57) 

From the boundary conditions fl7r|) we obtain the quadrature fiuctuations of the refiected 





2ai„^7m2^2 

7?+Vi 

-2ain72-y^7m2 
72 2 +^"2 
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fields, according to 



6pr = - ^^"^^^ + V^^Pi (58) 
J-^^^^^ + (59) 



out 



Spin + Spy r- 

+ \/^lra2C>P2 (60) 

+ 



= --^^^^^^ + (61) 

Finally, at the output port of the beam splitter, field quadratures are recombined according 
to 

5pBs = ^(-5pr* + 5pr*cos^ + 5gr*sin^) (62) 

Sqss = ^(-5gr* + '^C*cos^ + 5p°"*sin^) 
V 2 

where 9 is the phase shift introduced by the interferometer unbalance. 



A. Symmetric interferometer 

We will solve the above equations in the case of identical cavities and balanced Michelson 
interferometer. Therefore we replace ■jmi = 7m2 = 7m, In = In = 7«, ti = T2 = r, and we 
set ^ = 0. 

The field quadratures fluctuations at the output port of the beam splitter can be expressed 
by the product between input fluctuation Xj„ and two vectors of coefficients PBsl^Xi] and 
(5_Bs[5xj], where 6xi is the generic input fluctuation, with i =1 to 10: SpBs = Pbs " ^in and 
Sqbs = Qbs ■ Xin. The vectors Pbs and Qbsi obtained from Eqs. (15511621) . are the following: 

1 r ... /I , 0-4 



A ■ PBs[5pin] = — -n - Q4:k -fh ai„7„ {xiX2 - Xi2X2i)V'iV'2(V'i - ^2) 

(72 + ^0 (72 + L 

+ Seml^J [xi^i(f + - irn)i,l + 7((7 - ^rn){2^ - ztQ) + ^2')] 

- X2V'2(7' + i^lfbi^l + (7 - ^ra){2J^ - zjtQ + ^Pl)] 
+ S^hal-fJij^ + V>2)(72 + ^2)^27 - trQ)[xi2Ml{l " ^rQ) + ^l) 

- X2iV^2(7(7 - i-r^) +^?)] 
+ 7m(7 - ^rml' + V^?)'(f + i'lfii'l - 
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A-Pbs[Sp 



V 



X 



X { - 6Ak'h'ai^JixiX2 - Xi2X2i)V5iVS2[(7' + i^Dil' + 4') 

+ Sk'hal^^ [X2M^' + V^i)'[-77m((7 - ^rnf + ^l) + (-(7 - tr^f 

+ (7 - irnhm - + ^l)] - xiMi' + + mil - irnf + i^l) 

- lm{iir^{l - irVt) + (7 - irn)V'i - 7^2)] 

+ Sk'hal'jJi'y' + + m^l - trQ)[x2iMlil ' ^rQ) + i^f) 

+ xi2Mi{i - irn) + ^i)] 

- + mil' + m [((7 - ^rnf + mil - + 

- 7^(7 - iTO)[2(7 - IT^f + + i^l^ } 



A • PBs[5qv. 



ii'+mii'+^if 



X 



X 



I - 6Ak'^h'^af^-frn^ {X1X2 - Xl2X2l)lpl1p2ilpl - '02)(V'lV'2 - 7^) 



- 8k'hal-fJ [xiMi' + i'lfb'i'2 + m2 - mi - + i,l 

- X2M1' + ^D'b'^i + i^ii^l - mi - irnf + i^l)] 

- m''nmJmxi2 - X2i)^i^2(7' + mi' + ^2)(27 - 
+ 7m(f + mil' + i'lfii'i - mil - - m2)] 



A • PBs[5qv] 



ii'+mii'+^if 



X 



X 



iXlX2 - Xl2X2l)lpl1p2ilpl + '02)(7^ + 



- Sk'^haf^'jm 



xiMi' + mb'^2 + m2 + mi - ^m + ^d] 

+ X2M1' + ^'ifb'^i + i^ii^l + mi - irnf + v^Dlj 

- 8iehal^m''rn{xi2 + X2i)^i^2(7' + i^Dil' + ^2)(27 - irQ) 

- imii' + v^2)'(7' + mi^i + mil - ^m + m2)] 
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Pbs[Spii] 



(72 + ^2) (^2 + ^2) L 

+ 9>ehal^-fmMil{l'' + V52)((7 - irnf + V^') - X2V52(7' + i^lfil - ir 

- 8k''hal^rnX2iMl' + V5i)(27 - ^rfi)(7(7 - tril) + 

- (f + V^2)(f + V'i)'((7 - ^rnf + (7 - ^r^))V^2')} 



'12 



-Pbs[5pii] (1 ^ 2) 



A • PBs[5qii] 



Vl - 1mV^7^i>l{T + ^2 



(^2+^2)2(^2 + ^2 



If 
2 , J.2\ 



-|64A;Vaf„7„^(XiX2 - Xi2X2i)V'i'02 



+ Senal^^lxiMl' + V^2)((7 - + + ^^^^^^^ + ^2^2^ 

+ 8ik^halTnjrnX2iMl^ + V'i)(27 - ^7-0) 
+ (f + VS2')(f + VSD'((7 - irn)' + ^1)} 

PBs[Sqi2] = -PbsHi] (1 ^ 2) 



(^2 + ^2) (^2 + ^2 



8/c^^q;-„7„j(X2 + X2i)V'iV'2(27 - irQ) 



+ V'i(27 - iTn){^^ + V^2)((7 - ^rnf + 4')} 
^^55^02] = -PbsNoi] (1 ^ 2) 



(^2 + ^2)2(^2 + ^2)2 



|64A;S^^^ai„7m^(XiX2 - X12X21) 



+ Sk'hal^J [x2[7'(7' + ^i)'((7 - + 

- VSiV;2(f + V^D'lf + i^l)] - Xi[f (f + V'2')'((7 - irnf + i^l) 

+ 8A:2/,,«2 ^7^(72 + ^2)(^2 + Y;2)[x2i(7(7 " ir^) + i'Dilil - ^rQ) - ^ 

- Xi2(7(7 - ir^) - V^i)(7(7 - irQ) + i^i)] 

+ 7m(7' + + i^lfii^i - - irnf + } 
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+ Xi[7'(7' + ^2)'((7 - ^rl])' + ^^2) + ^1^2(7' + ^?)(7' + ^if] 

- 8k'hal^J[x2i{l' + i^Dil^ + ^2')(7(7 - ir^) + i^Dilil - i^) - i^l) 

+ Xi2(7' + i^Dil^ + i^Dilil - irn) - V;?)(7(7 - irO) + i^l)] 

+ 7m(7' + ^i)'(7' + + ^2)((7 - ^rnf + V^i^2)} 



2, , , -,,2{64A;S^^Q;f„7m^(xiX2 -Xi2X2i)V'i'02(V'? -"02) 



(^2 + ^2)^(^2 + ^2) 



+ Sk-'hal^J XiMl" + 4')'[(7 - iTQ)(i7TJl + - li^l] 



+ X2^2(7' + i^rbi^l ^ (7 - irn){i^Tn + ^2')] 

+ Stk'rnhal^Jh' + ^Dil' + ^i)[Xi2V^2(7(7 - irn) - ^l) 



A-QBs[Sqv] = 



(^2 + ^2)2^^2 + ^2 



X 



X ^-QAk n a,„7„ (xiX2 - Xi2X2i)^i^2 

X [(7' + V'i)(7' + ^2) - 77m(27' + + i^l)] 

+ Sk'hal^m [ - X2M1' + + i^Diil - irnf + i^l) 

- Imbi'l + (7 - irn){2-f^ - iirn + i^l))] 

- xiMi' + i^ifih' + ^i)((7 - ^Tnf + 

- iMl + (7 - ^TQ)(272 - i^rQ + i^D) - irQii^l - i^l)] 
+ SternUal^Ji^^ + ^^)(^2 ^ ^2)^^^^^^(^(^ _ ^^^) _ ^2) 

+ X2iV'i(7(7-«'r^^) -'i/'D] 

- ^^(7-^r^])(2(7-^r^])' + ^^^^2')]} 
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A n \x ^ V7-7mV^7mi7 + V2 ( aAiA 2*2-4 2/ w 

A • QbsWh] = — 7-27^^ - 64A; 7 ^ "m7m (X1X2 - Xl2X2l)^2 

(72 + -0f ) (72 + ^2) L 

.2fe-2 r , _.2/.,2 , „7;2\//., • r.N2 , „7;2\ , J. J. /.,2 , T2\2 



- ^k'haiimbcirir + V'2')((7 - ir^Y + i^i) + x2^^iV'2(y + v,r] 

+ 8ehaljmX2ih^ + ^i)(7(7 - ^^f^) - V^2)(7(7 - ^rQ) + 
+ ^^1(7' + V'2)(7' + ^l)'((^7 + mf - i^D] 

Qbs[5pi2] = -Qbs[Spii] (1 ^ 2) 



A • QnsiSqn] = ^7 , ^S'r^ J2I I " 64^V^«t7.^(XiX2 - Xi2X2i)V^i^2 

{i^ + Vi) (7^ + ^2) ^ 

- Senal^„^Mll{l^ + ^2')((7 - ^r^f + ^2) + X2^2(7' + - ^r^l)] 

- + + ^ifiil - irn)' + (7 - irn)i^l] } 

QBs[5qi2\ = -QBs[Sqii] (1 ^ 2) 

(7^ + A)(7^ + ^2) 



X 



(f + r2){{l - irnf + V;2')(7(7 - ir^) - i^l)} 



Qbs[5M = -Qbs[5M (1 ^ 2) 
where A is the determinant of N, expressed by 

^ ^ (^2 + ^2)^(^2 + ^2) {64fe%^«f„7^(XlX2 - Xl2X2l)^lV^2 

+ (7' + V^i)[(7 - irQ)^ + + Vi2')[(7 - irQ)^ + i^l] + 

8^2;ia^„7m [X2VS2(7' + VS?)[(7 - ^tQ)' + i^l] + 
Xi^i(7' + ^2)[(7-^rl])V^^j 



The homodyne detection after the beam sphtter allows to choose the output quadrature 

5Eout = ^Pbs cos w + 6qBs sin w, according to the detection phase w. It is useful to describe 
also SEout as product Vout ■ ^in, with a coefficients vector Vout defined as 

Vout = Pbs cos w + Qbs sin w . (63) 
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The input fluctuations ^■^oi and 5ijjo2 contain the gw signal. They also include ther- 
mal noise and any kind of classical fluctuations of the cavities length, that we are not 
considering in this article since in our approach they are not distinguishable from the 
gw signal. Taking an optimal interferometer orientation, the gw signal is proportional to 
d{Li — L2) = {2k)~^{dipoi — S1P02). Therefore we define a sensitivity Silfl) as 



(64) 



{2k)^\Vout[Si^oi]-Vout[S^o2]\' 

where Sx^xi are the spectral densities of the input fluctuations 5xi, which are assumed uncor- 
related. We take all noise spectra doublc-sidcd and normalized to shot noise. Therefore, for 
all the vacuum fluctuations (i.e., for 5pv, Sqv, Spn, Sqn, 6pi2, Sqi2) the spectra are SxiXi = 1- 
In totally symmetric conditions, we have ■0i = "02 = "0, Xii = X22 = Xs and Xi2 — 
X21 — Xc- To write clearer expressions, we deflne ^ = ■0/7 (detuning normalized to the half 
cavity linewidth); Pj„ = hkcal^ (input laser power); = Im/l] ^cav — l/^ (cutoff angular 
frequency of the cavity). The expression of Sl{^) is 

-1 



^2 



(2^'cosw + (*2 - 1) sinw) + 



n 



Xc 



3^2 _i ^(^2_3) 

1 H —r- cos zw H — — sm 2w + 



^ [i?e(xa-Xc)(*(l + ^') + ( 



sin w — ^ cos w) 
^2 



(65) 



+ 



+ 2^r„ 



2r. 



2^2 _ 2 - 



cos2w 



+ r„ 1 + - *2 6 + 



V ^cav J I V ^cav J I ' 



+ /m(x.-Xc)(2*(l + *')^(r^-l)) 



^4 ^ 



+ 1 + 



We notice that, in this case, Kut[(^'0oi] = Kut[(^'0o2] and the detection is not sensitive to 
the common mode of the interferometer. In the short cavity regime, with n/ilca^; — 0, the 
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Eq. (165|) reduces to 



+ -^Reixs - Xc) k(l + + r„ sin 2w (l + - 6^') 




Tm P [2^ COS w+ - I) sin w] ^ 




2r 



^2 - 



(66) 



+ ATm'i^i^^ -l)cos2w +(1 + ^2)^ 
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